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ABSTRACT 


This  report  contains  a  summary  of  it  .«  activities  carried  out 
under  Signal  Corps  Contract  No.  DAA.B07-e.H-c-  •  38 1  (Texas  A  i  M 
Research  Foundation  Project  No.  5R6)  during  rhe  cnnf.icr  period 
15  December  I960  to  15  June  1970. 

Activities  during  this  period  hove  been  devoted  primarily  to 
the  programming  of  the  equations  for  the  forested  boundary  layer 
and  the  wiring  and  checking  of  rhe  computer  patchboards  for  solu¬ 
tion  of  these  equations  on  the  general  purpose  analog  computer  at 
Texas  AMI  University.  The  complete  set  of  these  equations  scaled  in 
analog  format  is  included  in  this  report  as  are  wiring  diagrams  for 
solution  of  the  scaled  equations. 
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I. 


Introduction 


An  experimental  set  of  equations  for  simulating  the  atmosphere 
within  and  above  a  tropical  forest  is  under  development  at  Texas  A&M 
University.  The  general  analog  format  of  these  equations  is  presented 
in  Appendix  B  of  Technical  Report  F.CC+f-0381-3,  April  1970,  under 
Contract  DAA507-68-C-0381 .  The  complete  set  of  tentatively  scaled 
equations  for  the  tropical  forest  system  was  set  forth  in  that 
report.  The  time  scaling  used  in  the  initial  simulation  was  a 
computer-time  versus  real-time  factor  of  1940  to  1;  that  is,  the 
computer  solution  is  1490  times  faster  than  real  time,  so  that  a 
12-hr  solution  is  obtained  in  30  sec. 

The  principal  activities  during  this  reporting  period  have  been 
devoted  to  the  task  of  modifying  some  of  the  scaling  factors  initially 
employed  in  programming  the  equations  for  solution  on  the  General 
Purpose  Analog  Computer  at  Texas  A&M  and  to  wiring  of  the  computer 
patchboards. 

Section  II  of  this  report  contains  the  modified  scale  factors 
and  constants  employed  in  the  forest  and  boundary  layer  simulation. 
Section  III  contains  a  complete  set  of  equations  scaled  for  the  nev 
scale  factors,  and  Section  TV  contains  the  analog  computer  program, 
for  solution  of  the  equations. 
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IT.  Scale  Factors  and  Constants  Lmployed  in  the  Forest  Boundary 
Layer  Simulation 

A  complete  list  of  scale  factors  and  constants  employed  in  the 
forest  boundary  simulation  *s  shown  in  Table  1.  A  significant 
modification  has  been  made  to  the  scale  factors  for  the  east-west 
and  north-south  components  of  the  wind  in  the  free  air  above  the 
forest.  These  wind  components  were  initially  scaled  for  50  m/sec 
in  the  free  air  and  10  m/sec  within  the  forest.  This  value  of  50  n/sec 
is  certainly  far  in  excess  of  what  would  normally  be  expected  in 
cropiccl  areas,  except  in  very  extreme  circumstances,  such  as  in  a 
hurricane  or  typhoon;  therefore,  the  east-west  and  north-south  com¬ 
ponents  of  the  wind  and  free  air  have  been  rescaled  to  25  m/sec.  The 

> 

value  of  10  m/sec  within  the  forest  has  been  retained.  For  dense 

tropical  forests,  this  value  is  considered  to  be  extremely  hiph; 

however,  for  relatively  open  forests,  the  value  is  not  considered  to 

» 

be  excessively  large,  especially  near  the  edges  of  such  forests. 

The  east-west  and  north-south  components  of  the  shearing  stress 

2  2 

have  been  rescaled  similarly  from  20  dynes/cm  to  10  dynes/cm  .  This 
value  is  more  consistent  with  the  rescaled  value  for  the  east-west 
and  north. -south  components  of  the  wind.  The  other  scale  factors 
remain  unchanged,  as  do  the  fixed  constants  shown  in  Section  B  of 
Table  1.  These  scale  factors  also  can  be  considered  as  tentative 
since  modification  may  be  anticipated  subsequent  to  obtaining  an 
initial  solution  of  the  equations  on  the  general  purpose  analog 


computer . 


Table  1.  Scale  f^otors  and  constants  employed  in  flu*  forest 
boundary  layer  simulation. 


A.  _ ScaJ_e 

Factors 

I\.  : 

in  cm/sec 

R  : 

*• 

n 

.05  ii/cr  ^ec 

e  : 

100  mb 

s,  * 

cm/sec 

1  ,F  ; 
x  y 

10  dynes/cm'* 

SA.  : 

l-uor  cr/vec: 

FCC  : 

100  cm /km 

T  : 

I  a.  : 

100  cm/kn 

time  : 

lA4f>:l 

«ft 

n 

K  : 

r 

50.000  cm* /sec 

u,v  : 

1’ 500  cm/sec  (free  air) 

\  : 

.05  cal/cm' sec 

’'•00  cm/sec  (forest) 

c  : 

e 

2 

.05  cal/crr.^sec 

V  V 

2 

10  dynes /cm 

B.  Tons t ants 


a 

ar 

.621 

r  = 

K!  ,4t'  1  sec 

C 

— 

.24  cal/gn  deg 

p  = 

1012  mb 

P 

R 

= 

'J 

980  cr./sec" 

7.  ~ 

105-b  r  j 

I 

s 

*> 

.03)0  cal/cm.~sec 

j  - 

1.354x10  ^  cal/cm 

k 

K 

.40 

u>  * 

7.3x10  ^  rad/sec  ■ 

L 

. 

503.3  cal/gm 

p  * 

1.2x10  ^  gm/cm 
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III.  Scaled  Analog  Format  of  the  Equations  Employed  in  the  Forest 
Simulation 

v 

Table  2  contains  the  complete  expanded  set  of  the  newly  scaled 
equations.  Colu.nns  1  and  2  in  Table  2,  headed  j  and  1,  are  indices 
which  refer  to  the  particular  layer  of  applicability  of  the  equation 
involved.  The  third  column,  headed  z.  indicates  the  nominal  height 
in  meters  of  the  layer  to  which  the  individual  equation  applies. 
Similarly,  column  4,  headed  z^,  confnins  the  applicable  height  in 

r  , 

meters  of  the  layer  interfaces  at  which  the  fluxes  of  momentum,  heat, 
and  vapor  pressure  are  calculated,  and  column  5  contains  the  scaled 
equations  themselves.  The  numbers  in  parentheses  appearing  on  the 
extreme  right  of  the  equations  are  reference  numbers,  corresponding 
to  the  equation  numbers  given  in  Semi-Annual  Report  No.  3  on  Project 
586,  Contract  No.  DAAB07-68-0381. 

^  The  complete  set  of  equations  is  shown  here  to  facilitate  under¬ 
standing  of  the  computer  solution  diagrams  shown  in  Section  IV.  The 
coefficients,  c,  for  these  equations  are  tabulated  in  Table  3  for 
easy  reference.  Seme  of  the  values  of  c  appearing  in  this  table 
differ  from  those  given  in  the  aforementioned  report.  These  differences 
have  resulted  from  the  selection  of  new  scale  factors  for  the  east-west 
au%  north-south  components  of  the  wind  and  for  the  shearing  stresses. 

The  scaled  equations  have  been  written  so  that  the  scale  factors  are 

shown  explicitly  for  each  variable.  For  example,  in  equation  (si), 

\ 

at  the  level  of  250  m,  the  east-west  component  of  the  wind  is  scaled 
for  a  maximum  value  of  2500  cm/sec  and  the  time  scale  factor  is 


1440:1;  thus,  the  equation  is  divided  by  2500  4^d  multiplied  by  the 
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Table  2,  Continued 
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Table  3.  Korm  and  numerical  values  of  the  coefficients,  c,  in  the  scaled  equations 


Table  3.  Continued 
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time  scale  factor  of  1440.  The  scale  factor  for  A(g3h/3x)  la  shown 
to  be  1.736.  This  is  the  wind  scale  factor  of  2500  divided  by  a 
tine  scale  factor  of  1440.  The  advection  tens  (a^u^  -  87*7) 
are  shown  to  be  scaled  over  1.736.  The  first  ten  on  the  right  hand 
side  of  the  equation  is  written  as  .4g(3h/3x)/.6944.  Th<*  ten  is 
written  in  this  Banner  to  show  the  one  to  one  correspondence  between 
the  equation  fomat  and  the  calculation  of  this  ten  in  the  analog 
computer  solution  diagraau.  In  the  solution  diagrams,.  g(3h/3x)g/. 6944 
is  calculated  and  then  nultlplied  by  .4  to  obtain  a  scale  value  for 
the  pressure  gradient  ten  of  1.736. 

The  north-south  component  of  wind  in  the  Coriolis  ten  is  scaled 
for  a  maxi bub  value  of  2500  cn/sec.  The  coefficient  c^  is  shown  in 
Table  3  to  be  1440  (2u)  sin  #  and  applies  at  each  atmospheric  height. 

In  addition,  the  east-west  component  of  the  shearing  stress  of  the  wind 

2 

is  scaled  for  a  maximum  of  10  dynes/cm  .  The  coefficient  c*  _  is  shown 
in  Table  3  to  consist  of  the  time  scale  factor  divided  by  the  wind 
scale  factor,  multiplied  by  the  scale  factor  for  the  shearing  stress 
divided  by  air  density  and  the  thickness  of  the  layer  centered  at  250  m, 
and  further  multiplied  by  the  ratio  of  the  momentum  exchange  coefficient 
at  300  m  divided  by  the  momentum  exchange  coefficient  at  200  m  minus 
one.  The  calculated  values  for  each  level,  j,  is  shown  in  Table  3. 

The  equations  in  the  set  which  are  not  solved  explicitly  on  the 
Ceneral  Purpose  Analog  Computer  are  listed  in  unsealed  format  in 
Table  4.  These  equations  consist  of  equation  (5)  which  describes  the 
lapse  of  pressure  with  increasing  height;  equation  (11)  which  relates 


the  specific  humidity  to  the  vapor  pressure  and  the  atmospheric 
pressure;  equation  (12)  which  expresses  the  rate  of  gain  or  1 ot.s  of 
heat  per  unit  mass  to  the  system;  equation  (13)  the  equation  of  state 
for  a  perfect  gas;  and  equation  (14)  which  relates  the  vertical 
exchange  coefficients  for  momentum  heat  and  water  vapor.  Since 
these  equations  appear  Implicitly  in  the  wiring  diagrams,  scaling  is 
not  required. 

Two  maabers,  one  of  which  is  followed  bv  the  letter  f  (forest), 
are  assigned  to  each  equation.  This  double  numbering  has  been  done 
in  order  to  point  out  the  fact  that  the  equations  apply  to  both  the 
free  atmosphere  above  the  forest  and  the  atmosphere  within  the 
forest. 

Table  4.  Implicit  relations  in  the  solution  diagrams. 


PQ  -  Pgz 
ae/p 

C  dt 
P 

p/R  T 

a 

\-Kv 


(5)  and  (5f) 

(11)  and  (130 

(12)  and  (I4f) 

(13)  and  ( 1 5 f ) 


(14)  and  (16f ) 


IV.  Anal oc,  Solution  Diagrams  for  Tropical  Forest  Simulation 

? 

i  * 

There  are  many  ways  to  solve  the  scaled  equations  on  the  analog 
computer.  The  methods  employed  depend  upon  many  factors,  but  are 
primarily  controlled  by  the  quantity,  quality,  and  nature  of  the 
available  %6mputing  components.  The  method  shown  here  is  only  one  of 
many  that  could  be  employed  and  would  be  compatible  with  the  available 
computing  equipment. 

Efforts  have  been  made  in  setting  dp  the  solution  diagrams  to 

\ 

* 

treat  each  parameter  as  a  variable  as  far  as  is  possible.  The 
valuq/1  of  these  variables  may  be  easily  changed  from  case  to  case, 
although  for  a  particular  case,  the  variable  may  be  held  constant 
during  the  solution  period.  Discussion  of  the  diagrams  on  a  one-to- 
one  correspondence  basis  with  the'  equations  shown  in  Section  III  is  not 
practical  since  there  is  considerable  duplication  of  the  wiring 
patterns  due  to  the  similarity  of  the  equations-  for  the  ‘individual 
layers;  however,  a  brief  explanation  of  the  solution  of  each  type  of 
equation  will  assist  the  reader  in  making  the  correspondence  for  the 
solution  of  the  individual  equation?  for  the  various  layers. 

The  symbols  employed  in  the  analog  diagrams  and  their  inter¬ 
pretation  are  shown  in  Table  5.  The  left-hand  column  of  this  table 
shows  the  symbol  employed  and  the  right-ham’  column  indicates  its 
interpretation. 

The  analog  solution  diagrams  begin  on  page  30  on  which  the 
computer  diagrams  for  the  solution  of  the  east-west  wind  components 
appear.  This  diagram  is  the  programmed  solution  of  equation  (si) 


Table  5.  Analog  wiring  symbols. 


High  gain  amplifier 
with  an  input  to  the 
summing  junction,  an 
input  resistor  of  .1 
megohm,  an  input 
directly  to  the  grid. 


and  no  feedback 
Inverter 


B 

SIN 

-B 

COS 

Servo  resolver 


cup 


Capacitor 


High  gain  amplifier 
with  one  r.-vgohr  feed¬ 
back  resistor  and 
inputs  of  1 ,C  and  .1 
megohm  (surcer) 


Variable  function 
generator 


/"'N _ Potentiometer  (low 

side  grounded) 

Iiode 


/ 


multiplier 


rk 


f  round 


Follow-up  cup  of  a 
servo  multiplier  or 
resolver 


Function  switch 


H 

L 

CT 

T 


Servo  multiplier 
cup  (other  than 
fol lou-up) 


Low  speed  re  lav 
comparator 


RIAL 
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EAST  -  WEST  COMPONENTS  OF  WINO' 


shown  on  page  5.  The  solution  for  the  east-west  component  of  the  wind 
for  the  layer  at  150-m  height  appears  in  the  lower  left-hand  corner  of 
page  30.  ’ 

i 

The  east-west  wind  component  -u^q/2500  is  obtained  as  the  output 

of  an  integrator  which  sums  the  individual  terms  appearing  on  the  right- 

hand  side  of  equation  (si).  The  initial  condition  of  the  east-west 

wind  component  at  150  m  is  obtained  from  observed  values  and  is  sup- 

* 

plied  to  the  Integrator  through  the  initial  condition  potentiometer. 

This  potentiometer  supplies  the  initial  voltage  from  which  the  ampli¬ 
fier  integrates.  The  remaining  inputs  are  calculated  in  the  following 
manner. 

1.  Vertical  divergence  of  the  shearing  stress.  The  vertical 

divergence  of  the  shearing  stress  is  calculated  from  the  winds  at  75, 

150,  and  250  m.  The  solution  diagram  is  shown  in  the  upper  center 

portion  of  page  32.  Voltages  representing  Hi^q/2500  and  -u^.q/2500 

*  *». 

are  applied  individually  to  10-gain  inputs  of  an  amplifier  having  a 

potentiometer  set  to  (Z_crt  -  Z, ,_)  K  /c/1.25  x  10^  p  K  -cn  as  a  feed- 

250  150  m,45  m,250 

back  element  to  calculate  -(t  Oft-/10)/(K  /c/50,000).  The  output  of 

Xy  ^UU  in,  4  j 

this  amplifier  is  multiplied  by  ^,./50,000  by  means  of  a  servo- 

multiplier.  The  product  of  these  terms  -t^  *8  obtained  from  the 

arm  of  the  servomultiplier  and  the  east-west  component  of  the  shearing 

stress  +t  ,_-./10  is  determined  in  an  analogous  manner  from  the  next 
x,  100 

lower  amplifier  and  servomultiplier.  The  terms  4-r^  100^^  anc* 

* 

All  amplifiers  employed  in  the  solution  diagrams  are  of  the  inverting 
type;  therefore,  the  terms  unity  gain  or  10  gain  imply  a  sign  change 
within  the  amplifier.  • 
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~Tx  200^^  are  a<^e<*  *n  a  su®n*n8  amplifier  to  obtain  the  difference 

of  the  stresses  at  the  two  levels  +  (t  ~nn  -  t  --_)/10.  This 

x,200  x,100 

difference  is  further  multiplied  by  5.76/p  AZ^q  on  a  potentiometer 
and  applied  at  a  gain  of  one  to  the  input  of  the  -u^q/2500  integrator. 

2.  Coriolis  term.  The  Coriolis  term  for  the  east-west  component 
of  the  wind  at  150  m  is  obtained  by  multiplying  the  north-south  wind 
component  at  150  m,  calculated  by  the  -v^^q/2500  integrator  on  page 

39  and  inverted  in  sign  by  a  unity  gain  amplifier,  by  2880  u>  sin  on 
a  potentiometer  and  applying  the  resultant  voltage  at  a  unity  gain  on 
the  -u^^q/?500  integrator. 

3.  Wind  advection.  Wind  advection  is  calculated  for  the  layers 
at  heights  of  1.5  m,  35  m,  45  m,  and  250  m.  Advection  between  these 
layers  is  obtained  by  linear  interpolation.  Advection  of  wind,  temper¬ 
ature,  and  vapor  pressure  is  handled  in  this  manner.  Since  linear 
interpolation  is  used  at  150  m  height  for  advection,  two  inputs  are 
required  to  the  -u^^q/2500  integrator*.  There  is  a  contribution  from 
the  advection  at  both  45  m  and  250  m,  the  calculation  of  which  appears 
on  page  34.  Positive  and  negative  values  of  the  east-west  and  north- 
south  wind  components  are  applied  independently  to  two  function 
switches.  The  up  position  of  the  first  function  switch  is  connected 

to  +U2^q/2500  and  the  down  position  of  the  switch  to  -u^^q/ISOO.  The 

up  position  of  the  second  function  switch  is  connected  to 

and  the  down  position  to  +v25q/ 2500.  The  arm  of  the  firsft  function 

switch  is  connected  to  a  potentiometer  having  a  setting  of  1440 

thus,  the  voltage  appearing  on  the  ana  of  this  potentiometer  represents 

+1440  <*250  u250^^^  °r  °’50  U250^**^  depending  upon  whether  the 
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function  switch  is  in  the  up  position  or  in  the  down  position. 
Similarly,  the  up  position  of  the  second  function  switch  calculates 


-1440  825qv250^^^  at  t*ie  arm  t*ie  connected  to  the 

arm  of  the  second  function  switch.  These  two  terms  are  added  in  a 

t  "v 

summing  amplifier  to  obtain  the  advection  at  250  m,  ±(a'u. 


250 


3'V250)/1*736* 


The  neutral  positions  of  the  function  switches  are  wired 
individually  to  the  arms  of  two  potentiometers.  The  initial  values 
of  the  winds  are  set  on  these  potentiometers  so  that  advection  may 
be  held  constant  at  the  initial  value  throughout  the  solution  period. 
Whether  advection  is  held  fixed  at  the  initial  value  or  allowed  to 
vary  with  wind  speed,  the  gradients  of  the  wind  a'  and  £’  are  held 
constant  throughout  the  solution  period. 

Advection  at  45  m  is  determined  in  precisely  the  same  manner  as 
that  at  250  m.  These  advection  terms  are  multiplied  by  appropriate 
factors  appearing  on  potentiometers  of  which  the  arms  are  connected 
to  unity  gain  inputs  on  the  -u^^/2500  integrator  to  provide  linear 
interpolation  between  45  m  and  256  m  height. 

4.  Surface  height-contour  gradient.  The  solution  diagram  for  the 
gradient  of  the  surface  height  contours  is  found  at  the  top  of  page  34 
for  the  east-west  height-contour  gradient.  This  solution  diagram  is  the 
schematic  representation  of  equation  (s23)  on  page  11.  A  potentiometer, 
having  the  conq;uter  reference  voltage  applied  to  its  high  potential 
contact,  is  set  to  the  value  g  ICC  sin  (ICC)/6944j ICC  sin  (ICC)'.  This 
factor  is  multiplied  by  ICC  sin  (ICC)  which  is  set  on  the  second 


i- 
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of  Che  series  potentiometers.  The  arm  of  chls  second  potentiometer 

has  as  an  output  the  product  -g  ICG  sin  (ICG)/6944.  The  sign  cf  this 

product  is  determined  by  the  ratio  ICG  sin  (ICG) /| ICG  sin  (ICG) |  which 

% 

is  included  in  the  first  potentiometer  setting. 

The  product  -g  FCG  sin  (FCG)  is  computed  by  the  two  series 
potentiometers  which  appear  just  below  those  for  the  initial  contour 
gradient.  These  two  values  are  added  by  a  summing  amplifier  and 
multiplied  by  1440/3800  FI,  where  FI  is  the  prediction  time  interval 
expressed  in  hours,  and  applied  to  either  a  unity  gain  or  10-galn 
input  of  an  integrator  depending  upon  the  setting  of  the  potentiometer 
which  supplies  the  integrator  input  voltage.  If  the  potentiometer 
is  set  for  a  solution  period  greater  than  or  equal  to  4  hr,  a  gain 
of  unity  is  selected  by  placing  the  function  switch  in  the  up  position 
For  solution  periods  of  less  than  4  hr,  a  gain  of  10  is  selected. 

The  output  of  the  integrator  is  the  expression  occurring  under 
the  intep.ral  sign  of  equation  (s23)  save  a  factor  of  .1000  which 
appears  on  the  potentiometer  following,  the  integrator;  thus,  the  arm 

of  this  potentiometer  carries  the  value  1440  g  {(TCG  sin  (FCG)  - 

- 

ICG  sin  (ICG) 1/6944  Atldt,  where  At  is  in  seconds.  TVo  additional 
series  potent ioroeters  calculate  -g  ICC  sin  (ICG)/6944.  These  two 
values  are  added  by  a  summing  amplifier  to  obtain  +g  36/3x1^/. 6944. 
This  solution  procedure  permits  linear  variation  of  the  surface 
contour  gradient  with  time.  The  potentiometer  set  to  .1000  is  used 
simply  to  allow  higher  input  voltages  to  be  applied  to  the  integrator 
for  long  solution  periods. 
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The  east-west  component  of  the  surface  contour  gradient  is 
inverted  in  sign  by  an  amplifier  and  further  applied  to  a  high  gain 
amplifier  having  a  feedback  potentiometer  set  to  3600  (2u  sin  <f>) ; 


therefore,  the  gain  of  this  amplifier  is  1/7200  tu  sin  <f>  and  the 
surface  contour  gradient  is  multiplied  by  this  value  to  obtain  the 
north-south  component  of  the  geostrophic  wind. 

The  surface  height-contour  gradient,  -g(3h/3x)^/ .6944,  is 
further  multiplied  by  .4000  set  on  a  potentiometer,  of  which  the 
arm  is  connected  to  a  unity  gain  input  of  the  -u^q/2500  integrator 
and  which  scales  the  surface  height-contour  gradient  to  the  proper 
value  of  1.736. 

5.  Variation  of  the  contour  gradient  with  height.  The  variation 

of  the  height-contour  gradients  with  altitude  depends  upon  the  gradients 

of  the  mean  virtual  temperature  within  the  layer  under  consideration. 

The  computer  program  appears  on  page  32  in  the  lower  right-hand  corner. 

The  voltage  on  a  potentiometer  set  to  .3923  (3T/3x)  (ST/Sx).^^  is 

applied  to  an  inverting  amplifier  of  which  the  output  is  connected  to 

potentiometers  having  settings  of  a(?T/3x),c  +  b(3T/3x) „rft.  The'coef- 

45  250 

« 

ficients  a  and  b  have,  respectively,  the  vr.lues  .1356  and  .0107  at  75-m 

height,  .3810  and  .1312  at  150-m  height,  and  .5000  and  .5000  at  250-m 

height.  These  values  yield  output  voltages  which  represent  A [g( 3h/3x)  ]/ 

1.736  where  j  is  the  index  of  the  layer  in  question.  The  surface  contour 

gradient  is  assumed  constant  with  height  up  through  -45  m;  therefore,  this 

term  is  not  applied  at  and  below  45  m. 

At  the  heights  of  1.5,  15,  and  35  m  an  additional  input  is  applied 

t»  the  integrators.  This  input  represents  the  drag  upon  the  air 

produced  by.  trees  and  other  foliage.  The  calculation  of  this  drag 
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term  is  shown  on  page  40  and  consists  of  the  squaring  of  the  wind 

component  from  the  output  of  an  amplifier  by  applying  it  .to  the 

* 

x  and  y  inputs  of  an  electronic  multiplier.  This  product  is  inverted 

in  sign  by  use  ct  an  inverting  amplifier  and  further  applied  to  a 

4 

gain  o^  10  and  also  multiplied  by  the  coefficient  [(7.2  x  10  ACp)-l] 

v 

which  is  set  on  a  potentiometer  the  arm  of  which  is  connected  to  a 
gain  of  10. 

The  circuit  shown  in  the  lower  right-hand  corner  of  page  34  is 
a  timing  circuit  which  automatically  places  the  computer  in  HOLD  mode 
at  the  termination  of  a  computational  cycle.  In  this  mode  the 
solutions  are  stored  until  they  can  be  printed  by  a  digital  printer 
or  placed  on  punch  cards.  As  one  can  see  from  this  diagram,  +100 
volts  is  connected  to  the  normally  closed  contact  of  a  relay.  This 
voltage  is  applied  to  the  arm  of  the  first  set  of  contacts  on  the 
relay  and  subsequently  to  the  negative  terminal  of  this  relay,  RLAL. 
The  normally  open  contact  is  connected  to  +48  volts  as  are  the  arms 
of  the  remaining  two  sets  of  relay  contacts.  The  output  of  an 
integrator  is  connected  to  the  positive  terminal  of  the  relay.  The 
input  to  the  integrator  is  a  potentiometer  having  a  setting  of 
1440/ 3600  FT,  where,  FI  is  the  solution  interval  in  hours  when  the 
output  of  the  integrator  ts  less  than  +100  volts,  which  is  the  voltage 
applied  to  the  negative  contact  of  the  relay,  the  relay  coil  remains 
uncn^Riaed  and  the  contacts  on  the  relay  remain  in  the  normally 
closed  position:  however,  when  the  output  voltage  of  the  integrator 
exceeds  +100  volts  the  re  lav  is  energired  and  the  contacts  are 

38 


♦ 


% 


1 


r 


39 


NORTH -SOUTH  COMPONENTS  OF  WIND 


000*  » 


.switched  to  the  normally  upen  position.  Under  these  conditions,  +48 
volts  is  connected  to  the  negative  side  of  the  relay  input.  This 
arrangement  provides  a  positive  voltage  differential  of  approximately 


+52  volts  between  the  positive  and  negative  relay  inputs.  This  large 

1  .  •  •  Aide  i-!»  the  re  I  av  in.ni.it.  This 

differential  in  voltage  assures  a  positive  locking  of  the  relay  in 

-  :  cage  <!  i  f  i  crent  i a  1  of  approximate  ! 

the  energised  position  and  prevents  relay  chatter.  Simultaneously , 

'  ■  ivc-  relay  inputs.  This  nr  go 

+48  volts,  is  applied  to  the  normally  open  contacts  of  the  other  two 

;  i  -  nr-  :  a  ;>osi  t: 5  ve  j«vi  ing  of  the  re )  av  in 

sets  of  relay  contacts.  This  +48  volts  is  transferred  to  the  NA 

•  nr  'vent  •.  re  lav  rhater.  S  i '"ulr.arieousl  v , 

bus  of  each  console  through  isolating  diodes  in  order  to  energize  the 

■  ;  u'i'al  !  v  :  .•  n  ’''h!  acfr:  iif  the  otie'-r  ( wo 

hold  relays  which  stop  computation  at  the  time  specified  by  the 

■ ,  •  •  I  .  ■  *  .  '  !  ,  »  »  r  .  ;  r  i  j  ■  I  .  >  1- 1 . , .  ’ 

setting  of  the  potentiometer  connected  to  the  input  of  the  timing 
integrator.  The  output  of  the  timing  integrator  is  applied  to  the 
input  of  a  summing  amplifier  having  a  gain  of  1/2.  Tills  1/2-gain 
amplifier  simply  provides  a  convenient  means  for  monitoring  the  time, 
period  for  computer  simulations. 

The  solution  of  the  set  of  equations  representing  the  energy 
balance  at  the  top  of  the  forest  is  shown  on  page  43  and  consists  of 
wiring  diagrams  for  equations  sl9  through  s37  excepting  equations 
s23,  s24,  and  s32  through  s34.  The  solution  diagram  for  the 
calculation  of  the  cosine  of  the  solar  zenith  angle,  equation  s22, 
is  shown  at  the  top  cf  page  43.  The  solution  diagram  consists  of  a 
sine-cosine  oscillator  with  initial  conditions  for  the  sine  and  cosine 
of  the  local  hour  angle,  H.  The  cosine  of  the  local  hour  angle  is 
taken  from  the  second  integrator  and  multiplied  by  cos  4>  on  a 
potentiometer.  This  product  is  further  multiplied  by  cos  6  which 
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appears  on  a  series  potentiometer  to  obtain  cos  <f>  cos  5  cos  [15(t  -  12)], 
che  second  term  in  equation  (s22).  The  first  term  of  this  equation  is 
calculated  with  two  series  potentiometers.  The  first  potentiometer  is 
set  equal  to  sin  and  the  second  potentiometer  is  set  to  the  value  of 
sin  6  to  obtain  the  product  sin  $  sin  5  at  the  arm  of  the  second 
potentiometer.  These  two  terms  are  summed  in  an  amplifier  which  has  a 
gain  of  unity  for  positive  inputs  and  a  gain  of  zero  for  negative 
inputs;  therefore,  the  output  of  this  amplifier  is  -cos  t  for  the 
hours  of  daylight  and  zero  at  night.  The  output  of  this  amplifier 
is  applied  to  a  potentiometer  set  to  -N  so  as  to  calculate  e  ^  cos 
This  product  is  further  multiplied  by  the  solar  distance  factor  which 
is  calculated  through  the  use  of  a  potentiometer  as  an  additional 
feedback  on  a  summing  amplifier.  Multiplication  by  and  I  are  accom¬ 
plished  by  a  coefficient  potentiometer  set  to  and  an  amplifier  with 
its  feedback  modified  by  a  potentiometer  set  to  [(.05/1) —1 ] . 

Inspection  of  the  equations  for  the  free  air  |canopy  reveals 
that  each  equation  is  implicit  in  the  temperature  at  the  top  of  the 

canopy  T  /50;  therefore,  a  value  of  T./50  must  be  found  that  satisfies 
'  d  d 

* 

the  energy  balance  equation  at  the  top  of  the  canopy.  Since  the 
effect  of  T^/50  is  algebraically  the  same  in  each  equation,  the 
solution  diagram  is  equivalent,  analog  wise,  to  having  a  negative 
feedback  around  a  high-gain  amplifier  for  which  the  inputs  are  the 
energy  balance  terms.  This  condition  may  be  expressed  as 

*The  subscript  d  is  equivalent  to  the  height  of  40  m  on  the  wiring 
schemat i cs . 

\ 
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where  K  is  ths  open  loop  gain  of  the  amplifier. 

If  there  is  a  state  of  balance  for  this  feedback  network,  there 

is  a  solution  for  this  set  of  equations  in  terms  of  T,.  The  high 

d 

gain  of  the  amplifier  will  force  the  sum  of  the  inputs  to  be  zero 

and  the  amplifier  output  will  be  -T^/50. 

There  are  three  feedback  loops  for  the  -T./5Q  amplifier. 

d 

1.  The  output  of  the  -T^/50  amplifier  is  applied  to  the  neutral 
position  of  the  function  switch  and  further  to  the  input  of  the 

4 

patchboard  function  generator  which  calculates  +oT./.015.  The  out- 

d 

put  of  this  function  generator  is  applied  to  a  potentiometer  having 

e/3.333  as  its  setting.  The  resulting  product,  +ioTj/.05,  is  added  to 

d 

4  _ 

-oT^(m  +  lOn/e^/lOO) /.05  and  -(1  -  J)Sj/.05  to  obtain  the  net 
radiation  at  height  d.  ^/.05  is  further  fed  back  into  the  -T^/50 
amplifier  to  complete  the  loop. 

2.  The  output  of  the  -T^/50  amplifier  is  applied  to  the  neutral 
position  of  a  function  switch  and  thence  into  the  input  of  a  variable 


diode  function  generator  which  has  as  output  +e,  /100.  The  saturation 

Q  t  $ 

vapor  pressure  is  fed  to  another  high-gain  amplifier  which  has  an 

additional  input  of  -q  ,  /100C.  This  additional  input  provides  the 

e ,  a+ 

polarity  constraint  for  saturated  or  unsaturated  conditions.  The 

output  of  the  high-gain  amplifier  is  -e^/100  which  is  further  applied 

i ' 

to  10-gain  inputs  of  two  high  gain  amplifiers  on  page  35,  which  have 

i 

as  outputs  +(q  /.05)/(K  .  /50.000)  and  -(q  ,  /.05)/(K  .  /50,0()0). 

e,o+  m,<o  e,d-  — ,4? 

These  values  are  multiplied  by  ^/50,000  by  use  of  servomultipliers. 
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The  resulting  values  of  +q  05  and  -q  .  /.05  are  inverted  with  a  1 

e ,  d+  e,d- 

s* 

unity-gain  amplifier,  and  both  values  of  the  vapor  flux  are  fed  back 

into  unity-gain  inputs  of  the  -T./50  amplifier. 

d 

3.  The  third  negative  feedback  loop  results  from  the  calculation 

of  the  convective  heat  flux.  This  feedback  loop  is  shown  on  pages  43 

and  51.  The  output  of  the  -T./5G  amplifier  is  applied  to  10-gain 

d 

inputs  on  two  hi^h-gain  amplifiers  appearing  on  page  51.  The  output 

of  each  of  these  amplifiers  is  fed  to  the  high  side  of  a  servomulti- 

plier,  the  outputs  of  which  are  +q  ,,/. 05  and  -q  ,  /.05.  These 

c , d+  c , d- 

two  values,  subsequently,  are  fed  back  to  the  Inputs  of  the  -T^/50 
amplifier  to  complete  the  feedback  loop. 

The  .005  mfd  capacitors  used  as  feedback  elements  around  the 
high-gain  amplifiers  which  compute  -T^/50  and  -^e^/100  are  of  no 
computational  significance.  These  feedback  elements  are  used  exten¬ 
sively  throughout  the  wiring  diagrams  for  noise  suppression.  The 
low  speed  relay  shown  on  the  right  of  the  diagram  on  page  41  is  used 

simplv  to  provide  feedback  on  the  -T  /50  and  -e  ,/100  high  gain 

a  a 

amplifiers  while  the  computer  is  in  the  balance  mode. 

The  circuit  diagram  for  calculating  the  momentum  exchange 
coefficients  at  3,  30,  and  45  m  height  is  shown  at  the  bottom  of 
page  43.  This  circuit  calculates  equations  s34  through  s37.  In  the 
lower  left-hand  comer  of  page  4  3  -u^^/2500  and  -v^^/2500  are  inde¬ 
pendently  connected  to  the  x  and  y  inputs  of  two  electronic  multipliers. 

2  6  2  6 

Thc'outputs  of  these  multipliers  -u^^/6.25  x  10  and  -v^/ 6.25  x  10 

2  6 

are  further  added  in  a  summing  amplifier  to  get  +S^/6.25  x  10  ,  the 
square  of  the  magnitude  of  the  .vector  wind.  The  square  root  c‘ 


2  6 

^6.25  x  10  is  obtained  by  an  electronic  multiplier  used  in 

square  root  mode.  The  output  -S^^/2500  of  this  electronic  multiplier 

is  applied  to  a  unity  gain  input  on  a  summing  amplifier  and  to  the  x  • 

input  of  another  electronic  multiplier.  The  quantity  +(S^,.  +  a)/2500 

is  calculated  by  adding  -S^^/2500  and  -a/2500  in  the  summing  amplifier. 

The  quantity  a  is  chosen  to  be  3  m/sec.  This  sum  is  further  applied 

to  both  the  x  and  y  inputs  of  an  electronic  multiplier  to  calculate 
2  -  6 

—  (S^,.  +  a)  /6. 25  x  10  .  This  value  then  is  inverted  in  sign  and 

applied  to  the  y  input  of  an  electronic  multiplier. 

The  temperatures  +T/c/50  and  -T  /50  are  added  In  a  summing 

amplifier  and  multiplied  by  Ag/3.125  x  10^9  set  on  a  potentiometer. 

The  resulting  product  is  applied  to  the  x  input  of  an  electronic 

2  6 

multiplier  in  divide  mode  which  has  as  its  y  input  -(S^  +  a)  / 6.25  x  10  . 
The  quotient  obtained  is  the  modified  Richardson  number.  The  modified 
Richardson  number  is  then  applied  to  the  input  of  the  function 
generator  containing  -6  /40  as  its  output.  This  quantity  is  applied 

IC  ,  4 

further  to  the  y  input  of  an  electronic  multiplier  having  as  its  x 
input  -S^^/2500.  The  multiplication  of  these  two  terms  yields  the 
momentum  exchange  coefficient  at  45  m  height. 

The  calculation  of  the  exchange  coefficient  for  momentum  at 
1.5-m  height  is  completely  analogous  to  that  at  45  ra.  Linear  inter¬ 
polation  is  used  between  1.5  m  and  45  m  to  calculate  the  values  for 


K  .  and  K  .. . 
m,3  ®, 30 

The  solutions  of  the  equations  for  the  forest  surface  section 


are  similar  to  those  for  the  air  canopv  section  except,  in  that  case. 
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the  equations  are  implicit  in  the  surface  temperature  of  the  ground, 
-T^/50.  This  progr^  appears  on  page  49. 

For  the  surface  section,  an  integral  exchange  coefficient  is 

employed.  The  diagram  for  the  calculation  of  this  parameter  is  shown 

on  the  lower  right-hand  side  of  page  49.  The  input  of  +R^/.25  is  fed 

to  the  function  generator  containing  p  /. 006  as  an  output.  This  output 

6 

is  further  applied  to  the  x  input  of  an  electronic  multiplier  having  as 
its  v  input  -S  J./1000;  the  resulting  product  is  -D  /10. 

1  •  3  1  •  !) 

A  simplified  soil  structure  is  employed  in  the  forest  surface 

section.  Its  analog  solution  diagram  is  found  on  the  left  edge  and 
near  the  center  on  page  49.  The  soil  surface  temperature  +T'/50  is 
obtained  from  an  integrator  having  inputs  of  the  soil  heat  flux, 

-q  / . 0 1 ,  multiplied  by  2>'iT/Pui  set  on  a  potentiometer,  +T’/50 

applied  in  a  negative  feedback  loop  to  a  potentiometer  for  which  the 
setting  is  2 ~/P,  and  the  mean  representative  temperature  for  the  soil 
also  multiplied  bv  2-/P  set  on  a  potentiometer. 

The  temperatures  are  calculated  from  equations  s3  and  s3f  and 
are  shown  in  diagrammatic  form  on  page  50.  The  solutions  are  completely 
analogous  to  those  for  obtaining  the  winds  discussed  earlier.  The 
radiation  term  R  is  calculated  in  the  lower  lelt  hand  comer  of  page 
51,  and  the  value  tR/. 03472  ;s  applied  to  potentiometers  feeding  the 
individual  integrators.  This  arrangement  permits  the  variation  of  R 
with  height  and  forest  structure.  Temperature  advection  terms  are 
handled  in  exact lv  the  sane  manner  as  are  the  wind  advection  terns. 
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A!R  TEMPERATURES 


FLUXES  OF  CONVECTIVE  HEAT 


ADVECTION  OF  TEMPERATURE 


The  diagrams  for  calculating  vapor  pressures  are  shown  on 
page  54  and  are  seen  to  be  completely  analogous  to  the  calculations 
of  temperature  and  wind  components.  The  only  term  that  is  signi¬ 
ficantly  different  from  the  others  is  the  moisture  term  which  is 
calculated  on  the  lower  left-hand  corner  of  page  55.  This  term 

I 

applies  only  within  the  forest  where  trees  and  brush  act  as  moisture 
sources . 

This  set  of  solution  diagrams  represents  only  an  initial  step 
toward  the  development  of  an  analog  system  for  simulating  the  behavior 
of  the  atmosphere  within  and  above  a  tropical  forest.  Equation  refine¬ 
ment  with  resulting  modifications  in  the  solution  diagrams  will  be 
required  since  little  factual  information  exists  concerning  many  of 
the  parameters  employed  to  describe  the  forest  and  functional  defini¬ 
tions  are  essentially  nonexistant. 

The  moisture  source  term,  M,  presently  being  employed  as  a  con¬ 
stant  requires  functional  definition  as  does  the  radiation  term  R. 

In  addition,  the  magnitude  of  the  drag  coefficient  is  not  known 
for  air  flow  below  the  forest  canopy;  therefore,  for  lack  of  better 
information,  a  value  is  being  employed  that  is  representative  of  flow 
over  rough  surfaces. 

The  exchange  coefficients  for  momentum,  heat,  and  vapor,  now 
defined  by  the  Deacon  wind  profile,  may  need  redefinition  for  forest 
application.  In  this  regard,  the  roughness  of  the  forest  canopy  will 
require  reformulation  in  order  to  be  compatible  with  Deacon's  wind 
profile . 
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VAF'O*  PatS: 


AOVECTlON-  OF  VAPOR  PRESSURE 


V.  Engineering  Activities 


During  the  past  six  months,  the  maior  engineering  accomplishments 
have  consisted  of  the  completion  of  the  new  test  jip  for  checking  the 
individual  computing  components  for  the  general  purpose  analog  computer 
^CPAC),  the  installation  of  a  completely  now  Grounding  svscer  for  the 
CTAC,  and  fhe  installation  of  a  new  enclosed  re  lav  rack  v’ich  will 
contain  the  five  diode  function  generators  present lv  in  use  and  t!  ree 
additional  ones.  Other  engineering  tasks,  not  as  large  hut  just  as 
important,  were  completed  during  this  report ing  period.  Two  trunk 
lines  were  run  from  Console  4  of  the  CPAC  into  the  electronics  shop 
for  use  with  the  new  test  jig.  All  of  the  integrators  on  the  CPAC 
were  adjusted  again  for  proper  rates  of  integration.  All  of  the 
servonultipliers  and  servoresol vers  were  checked  for  alignment  and 
for  worn  bearings,  gears,  and  potentiometers.  Three  of  these  servos 
required  rebuilding.  In  addition,  he  nccpssv-  relays  and  switches 
were  added  tc  the  CPAC  remote  control  box  to  prrrit  control  of  Console 
5  simultaneously  with  Consoles  1  through  4. 

Cork  was  begun  on  updating  and  drafting  of  wiring  diagrams  for 
project  equipment  in  use  and  on  setting  up  a  suitable  procedure  for 
indexing  all  wiring,  schematics  for  easy  filing  and  ready  reference. 

New  equipment  received  included  a  six  channel  Brush  recorder, 
a  BRPK  Teletype  Corporation  hirh  f-pecd  iapr  Punch,  a  fluke  model  Tib 
voltage  standard,  and  a  Fluke  mode!  t  V"V  digital  voltmeter. 
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angle  made  by  the  initial  contour  gradient 
vector  with  the  y-axis 

(deg) 

i.j.k 
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J 

albedo 

(nnn-d iocnsional) 

Kh 

exchange  coefficient  for  heat 

(cn//sec) 

K 

01 

exchange  coefficient  for  momentum 

(cm^/see) 

K 

V 

exchange  coefficient  t. water  vapor 

( ra^ /sec) 

k 
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(non-dimensional • 
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latent  heat  of  vaporization  of  water 

(cal  /gm't 

ln 

net  longwave  radiation 

t 

(cal /< m'sec) 

M 

aoisture  source  term  for  forest 
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empirical  radiation  factor 
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turbidity 
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n 

empirical  radiation  factor 
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atmospheric  pressure 
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energy  addition  per  unit  mass  from 
non-adiabatic  processes 
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qc 

convective  heat  flux,  positive  upward 

(cal / cm'sec) 
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convective  heat  flux  in  air  layer  just  above 
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(cal /cm*sec) 

qc.d- 

convective  heat  flux  in  canopy  just  below  the 
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-> 

(cal /re’sfc) 

% 

evaporative  heat  flux,  positive  upward 

(cal/cm‘sec) 
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a 
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(cal /cm* sec) 
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S 
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*• 
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o 
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T* 
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t 
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u 
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,v 


E 


x 


u 

?ero 


soil  boat  -lux,  positive  downward 
radiations’  coolinr  or  warming 
gas  constant  for  dry  air 
net  radiation 

nean  specific  humidity  for  a  subscripted 
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r cos  trophic  wind  components 

east-west  coordinate,  positive  eastward 
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l'S.nor,  top  of  boundary  laver 

vertical  coordinate,  positive  upward 

surface  roughness  length 

as  a  subscript  indicates  value  at  x  •  0 
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rradient  of  wind  alone  streaeline 


(cal /cr  sec) 

(deg/scc) 
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stability  r,irr>  .  r  Ci.c  .eacon  profile 

(r  jr-  d  mens  iona 1) 

ratio  of  wind  speed  to  radius  of  curvature 
of  streamline 

(H'f  S 

solar  declination 

(dec  > 

vertical  interval  between  man  1  t  ights  of 
ad intent  simulator  lavers 

( rr  ) 

incrorrrt  of  tir. 

(See ) 

nean  vertical  th  i  cV.  nc :  ot  tie  Inver 

i  ,-a ) 

emissivitv 

( r.( - c  imens  i  ora ! 

solar  zenith  angle 

(dep) 

Bean  potential  tenpertftifre  for  a  subscripted 
layer 

( deg  h  ) 

1/2  depth  of  air-canopy  layer 

(cn) 

depci.  of  forest  surface  layer 

(cm) 

voluaetric  heat  capacity  of  soil 

3 

(cal/cts  dec) 

therral  conductivity  of  soil 

(cal/cn  sec  deg) 

aoisturc  paraneter  for  the  canopT 

( c^ 1 / cb ^ s o c  rb) 

3.14.  .  .  . 

(non-d irensi vnal) 

1.2  x  10  \  air  density 

(pr./cm3) 

1.354  x  10  Stefan  Boltroann  constarc 

(cal/ce^sec  dec^) 

coaponent  of  t  in  x  direction 

(dynes/ cxs~  ) 

component  of  r  in  y  direction 

(dyncs/cc* ) 

latitude  % 

(dec) 

forest  transrissivi ty 

(r.cr-d  iccns'icnji) 

solar  distance  factor 

(ren-d iners i cnal ) 

7.3  x  1°  (angular  velocitv  of  earth’s 

rotation) 

( r jd/src ) 
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